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Abstract. Within the class of inflationary models, k-inflation represents the most general 
single field framework that can be associated with an effective quadratic action for the cur- 
vature perturbations and a varying speed of sound. The incoming flow of high-precision 
cosmological data, such as those from the Planck satellite and small scale Cosmic Microwave 

■ Background (CMB) experiments, calls for greater accuracy in the inflationary predictions. In 
this work, we calculate for the first time the next-to-next-to-leading order scalar and tensor 
primordial power spectra in k-inflation needed in order to obtain robust constraints on the 
inflationary theory. The method used is the uniform approximation together with a second 
order expansion in the Hubble and sound flow functions. Our result is checked in various 



limits in which it reduces to already known situations. 
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1 Introduction 

Inflation [1, 2] (for reviews, see Refs. [3-7]), which is currently the leading paradigm to 
describe the physical conditions that prevailed in the very early Universe, is now entering a 
new phase. With the advent of new high-accuracy cosmological data [8-21], among which 
are the Planck data [22], one can hope to obtain very tight constraints on the inflationary 
theory and even to pin-point the correct model of inflation. In order to achieve this ambitious 
goal, one must be able to compare the inflationary predictions to the data. The problem is 
that the inflationary landscape is very large [23] and that there is a whole zoo of different 
models making different predictions. Moreover, for many of these models, predictions can 
only be worked out by numerical methods. It is therefore not obvious how to extract model- 
independent constraints on the inflationary scenario. 

How then should we proceed? Clearly, one can approach the problem step by step 
and start with the simplest models. In other words, it seems reasonable to consider more 
complicated models only if the data force us to do so and tell us that the simplest models are 
not enough. Then comes the question of identifying these models. One can convincingly argue 
that slow-roll Single Field with a Minimal Kinetic term (SFMK) scenarios are the simplest 
inflationary models since they are just characterized by one function, the potential V(4>). In 
order to establish their observational consequences, a possible approach is to scan models one 
by one and calculate the predictions exactly [24-29], most of the time numerically [30, 31] 1 . 
This leads to an exact mapping of the inflationary landscape within this class of scenarios 
but, given that the number of SFMK models remains large, it would represent a huge effort. 
Another approach consists in developing a scheme of approximation allowing us to derive 
analytical, or semi-analytical, predictions. Although this is not always possible, such a 
method is available for the SFMK models and one can explicitly write a functional form 
for the primordial power spectrum of the cosmological perturbations [32], and even their 
higher order correlation functions [33-38]. 

In fact, one can enlarge the class of what we consider as the simplest models of inflation 
and assume that these ones are k-inflationary scenarios. K-inflation [39, 40] encompasses 
standard inflation and is more general since not only the potential but also the kinetic term 

1 See for instance http://theory.physics.unige.ch/~ringeval/fieldinf.html. 
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is now a free function. At the perturbation level, the action for the comoving curvature 
perturbation has a varying speed of sound and this describes all possible quadratic terms 
within the effective field theory formalism [41, 42]. But, more interestingly, and despite 
the fact that this class of scenarios is more complicated to analyze, a properly generalized 
slow-roll approximation can still be used. 

1.1 State-of-the-art 

At this stage, it is interesting to recall the present status of the techniques that enable us to 
calculate the two-point correlation function for the primordial cosmological perturbations. 

The spectrum of density perturbations during inflation was computed for the first time 
in Refs. [43, 44] and for the gravity waves in Ref. [45]. Then, in Ref. [46], it was realized 
that it can be evaluated exactly in the case of power-law inflation. The first calculation at 
first order in the so-called "horizon flow parameters" and using the slow-roll approximation 
was performed in Ref. [32]. This calculation was done for the SFMK models. This is a 
fundamental result since it allows to connect the deviations from scale invariance to the 
microphysics of inflation. This result was re-derived using the Green function methods in 
Ref. [47], using the Wentzel-Kramers-Brillouin (WKB) method in Ref. [48] and using the 
uniform approximation in Refs. [49, 50]. In fact, the Green function method of Ref. [47] made 
possible the first determination of the scalar power spectrum at second order in the "horizon 
flow parameters". Indeed, at second order, the mode equation describing the evolution of 
the cosmological perturbations can no longer be solved exactly, hence the need for a new 
method of approximation. Higher order corrections were also obtained in Ref. [51]. The first 
derivation of the tensor power spectrum at second order using the Green function method was 
presented in Ref. [52]. In Refs. [53, 54], it was also shown how to improve the WKB method 
by adding more adiabatic terms. This improved WKB method has allowed a re-derivation of 
the scalar and tensor power spectra at second order and confirmed the results of the Green 
function approach. 

After the advent of k-inflation, various attempts have been made to derive the corre- 
sponding power spectra. The problem is complicated due to the fact that density pertur- 
bations now propagate with a time-dependent speed (the speed of sound). In Ref. [55], the 
Green function method has been used but with some extra-assumptions on the behavior of 
the sound speed. The question was also considered in Refs. [56, 57] but the results obtained 
in those articles were not totally correct since the sound speed was (implicitly) assumed to 
be constant which is not the case in most of the k-inflationary scenarios (this result was also 
used afterward in Ref. [58]). These works also missed the influence of the sound speed in the 
tensor power spectrum due to the shift between the scalar and tensor pivot scales [59, 60]. 
The first fully consistent result for the k-inflationary scalar power spectrum was presented 
in Ref. [61]. The latter has been re-derived using the uniform approximation in Ref. [59] 
together with the first fully consistent calculation of the tensor power spectrum at the same 
pivot scale. These spectra were compared to Cosmic Microwave Background Anisotropy 
(CMB) data first in Ref. [62]. However, all of these calculations have been derived at first 
order only and no complete result at second order exists in the literature. 

The main purpose of this article is to close this gap and to derive the slow-roll power 
spectra for the density and tensor perturbations in k-inflation, at second order in the Hubble 
and sound flow functions 2 . This calculation is interesting for two reasons. Firstly, the second 
order result is available for SFMK models and, for completeness, it should also be done for the 

2 Conforming to the modern usage, we will prefer the denomination of "Hubble flow functions" and "sound 
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k-inflationary models. Secondly, according to the "blue book" [63], Planck will measure the 
spectral index with accuracy An s ~ 0.005. Even if one expects the Hubble flow parameters 
to be less than 10 -2 , second order corrections will be of order 10 -4 , that is to say relevant 
for high-accuracy measurements of n s and/or estimation of the corresponding error bars. 
Moreover, having at hand the second order terms allows to marginalize over, a procedure 
that should always be carried on to get robust Bayesian constraints on the first order terms. 

Before moving to the calculation, let us briefly recall some well-known results about 
k-inflation at the background and perturbation levels. 

1.2 K-inflation in brief 

K-inflation corresponds to a class of models where gravity is described by General Relativity 
and where the action for the innaton field is an arbitrary function, P((p,X), the quantity X 
being defined by X = —(l/2)g^ u d^(j)d u (j). This action can be written as 



b = — — I d xy/—g 



<P1 



(1.1) 



where M P1 is the reduced Planck mass. In fact, in order to satisfy the requirements that the 
Hamiltonian is bounded from below and that the equations of motion remain hyperbolic, the 
function P(X, <fi) must satisfy the following two conditions [64] 

dP d 2 P dP , , 

ax >0 - 2X l»r- + 3x >0 - (L2) 

The general action (1.1) includes standard inflation for which P = X — V ((/)), where V(<p) is 
the innaton potential. This class of model is in fact characterized by an arbitrary function of 
(f) only. K-inflation also includes the Dirac-Born-Infeld (DBI) class of inflationary models [65] . 
For those, one has P = — T((/>)y / l — 2X/T((p) + T(tfi) — V(4>). This kind of action typically 
appears in brane inflation and T{4>) is interpreted as a warping function representing the 
bulk geometry in which various branes can move. It is of course possible to find even more 
complicated examples but, in the following, we will not need to specify explicitly the function 

P{x,4>). 

As in standard inflation, the dynamics of the background space-time can be described 
by the Hubble flow functions e n defined by 

dlne n _ F ini 

e " +1 -^\T' e ° = lT' (L3) 

where N = \a.{a / a- ia \) is the number of e- folds. Inflation occurs if e± < 1 and the slow-roll 
approximation assumes that all these parameters are small during inflation e n -C 1. Let 
us notice that it is difficult to have an inflationary model without such a condition because 
otherwise one would obtain a deviation from scale invariance which would be too strong to 
be compatible with the cosmological data (see however Ref. [42]). 

At the perturbed level, we have density perturbations and gravity waves. As usual, rota- 
tional perturbations are unimportant since they quickly decay. Obviously, the tensorial sector 
of the theory is standard since the gravitational part of (1.1) is the ordinary Einstein-Hilbert 



flow functions" to refer to the original, but confusing, appellation "horizon flow parameters". See Sec. 1.2 for 
the definition of the Hubble and sound flow functions. 
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action. As a consequence, the equation of motion for the amplitude /ij, of gravity waves (re- 
scaled by a factor 1/a for convenience, where a is the Friedman-Lemaitre- Robertson- Walker 
scale factor) takes the usual form, namely 

ti£+[k 2 -U T (r,)]» k = 0, (1.4) 

where r\ is the conformal time and a prime denotes a derivative with respect to 77. The 
effective potential for the tensorial modes can be written as U T = a 2 H 2 (2 — ei), i.e. only 
depends on the first Hubble flow function (H = a' /a 2 is the Hubble parameter). 

For the density perturbations, the situation is slightly more complicated. One can 
show that the comoving curvature perturbation in Fourier space, can be written in terms 
of a modified Mukhanov- Sasaki variable vi~ by means of the following expression, v k = 
(ay / ii)Cfc/c s (in Planck units) where the quantity c s is defined by the following equation 

° l = Px + 2XP, XX ' (L5) 

a subscript ", X" denoting differentiation with respect to X. This quantity can be interpreted 
as the "sound speed" of density fluctuations. Notice that, because of the two consistency 
relations (1.2), we have c 2 > 0. The fact that c s is the sound speed can be most easily seen 
if one writes down the equation of motion of the Mukhanov- Sasaki variable. It reads 

v'^+[c 2 ( V )k 2 -U s ( V )]v k = 0. (1.6) 

This is similar to the equation of motion of a parametric oscillator. The quantity U s is the 
effective potential for the density perturbations and is a function of time only. As expected, 
Cg appears in front of the k 2 term, which is nothing but a gradient term in Fourier space 
and this confirms its interpretation as a time dependent sound speed. Since Cs(r]) is not 
known a priori, one can introduce a second hierarchy of flow functions in order to describe 
its behavior. Therefore, we define the sound flow functions <5 n 's by 

_ dln£ n _ c siui 

^ +1 = ^AT' 5o = — ■ (L7) 

Consistent models of inflation are obtained if 5 n <C 1, that is to say if the sound speed does not 
change too abruptly [59, 61]. A remark about terminology is in order at this point. In terms 
of the Hubble and sound flow functions, the effective potential for the density perturbations 
can be expressed as 



U s = a 2 H 2 



3 111 

ei + + ^ e 2 - 2 ei£ 2 + 2 €2€s + ( 3 ~ £ 1 + £ 2) #1 + S l + S l$2 



(1. 



The quantity U s depends on the e n 's up to £3 only and on the <5 n 's up to 82 only. Despite 
this last property, it is important to remember that the above expression of U s is exact and 
that no approximation has been made at this stage. 

The cosmological observables we are interested in are the two point correlation functions 
of the fluctuations, i.e. in Fourier space, the power spectra of both gravity waves and density 
perturbations: 



2k 3 

v h = — 

T* 1 



7T- 



2 Vt = — I4I 2 = — Cs2| " fc ' 2 (19) 
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They have to be evaluated at the end of inflation and on large scales. After the inflationary 
era, and for single field models, these power spectra remain constant and can be directly 
used to compute various observable quantities such as the CMB anisotropics or the matter 
power spectrum. Our goal is now to integrate the equations of motion of /i& and v \~ in order 
to explicitly evaluate the above power spectra. 

This article is organized as follows. In the next section, after having very quickly 
reviewed how the uniform approximation can be used in the cosmological context, we apply 
it to the calculation of the scalar and tensor primordial power spectra. Our results are 
discussed Sec. 3, in which we compare them, in the appropriate limits, with the existing 
literature and we present our conclusions. 

2 K-inflationary Power Spectra 
2.1 The Uniform Approximation 

In this section, we use the uniform approximation to calculate the power spectrum of the 
density fluctuations in k-inflation, at second order in the Hubble and sound flow functions. We 
have seen in the previous section that the density perturbations in k-inflation propagate with a 
time-dependent velocity Cs(ry). As the mode equation can no longer be solved exactly in terms 
of Bessel functions (even at first order for the sound flow functions), this prompts for the use 
of new techniques. Here, we choose to work with the well-suited uniform approximation [50]. 
The idea is to rewrite the effective potential according to U s = [y 2 — l/4:)/r] 2 , an equation 
which has to be understood as the definition of the function v(rj). Then, we introduce two 
new functions 



2/3 



(2.1) 



where the turning point time ?7*(fc) is defined by the condition g(j]*) = 0, that is to say r/* = 
—v(ri*)/[kc s (??*)]■ According to the uniform approximation, the Mukhanov-Sasaki variable 
can then be expressed as 

v k (v) = A k (?) 1 Ai (/) + B k (£) ' Bi (/) , (2.2) 

where the two constants A k and B k are fixed by the choice of the initial conditions and where 
Ai and Bi denotes the Airy function of the first and second kind respectively. Since one needs 
to compute v k on large scales, only the asymptotic behavior of the Airy functions is needed 
and one arrives at a simpler formula, namely 

j%^")=^^G /3/2 )' <2 - 3) 

Here, the function g(rf) should be taken in its asymptotic limit, i.e. g 1 / 2 ~ —v(j])/rj. Inserting 
the last equation for ^(77) into the formula (1.9), one obtains the following expression for 

k 3 \B k \ rjc 2 2is 

where we have defined \F = 2/ 3 / 2 /3. One verifies that is positive definite since the 
conformal time is negative during inflation. Therefore, the only thing which remains to be 
done is to express the combination c 2 /(a 2 uei) and the quantity at second order in the 
Hubble and sound flow functions. 
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2.2 Hubble and sound flow expansion 

The first step of the calculation consists in determining the functions a(rj), c s (ry), z/(r?) and 
€1(77) at second order in the Hubble and sound flow functions. Here, we first briefly explain 
the method in the case of the scale factor. By definition, the conformal time is given by 
77 = — J dt/a(t), where t is the cosmic time. By successive integrations by parts, one can 
re-write 77 as 



V 



1 

H 



1 + ei + e( + ei e 2 - aH 



1 d 

a da 



1 

H 



0? + e 1 e 2 ) 



da 



(2.5) 



where H = a' /a is the conformal Hubble parameter. It is important to stress that this equa- 
tion is exact. In the last term, the integrand is third order in the ej. Indeed, differentiating 
the term 1/H produces a e± which, multiplied with (ef + €±€2), is third order. We also have to 
differentiate expressions quadratic in the Hubble flow functions but, since de n /diV = e n e n+ i, 
this also gives third order quantities. Therefore, the last term is 0(e 3 ) and can be dropped 
for a second order calculation. In other words 



H 



[1 + ei + e? + eie 2 ) +0(e 3 ) 



(2.6) 



In fact, this equation is not exactly what we want yet because, although the second order 
terms e\ and eie 2 can be considered as constant in time 3 , this is not the case for the term 
ei which is first order. In order to render explicit the time-dependence, let us notice that 
the equations defining the Hubble-flow functions can also be written as de n /dr] = %e n e n j r \. 
Given that e n e n +i is already a second-order term, we can just replace % with — 1/?? in this 
expression and one gets e n = e n * — e n *e n+ i* In (77/77*) + 0(e 3 ), where we have chosen the 
integration constant such that this approximation is accurate around the time 77* of the 
turning point. Inserting this expression into Eq. (2.6) gives 



n 



- (1 + ei* + e?„ + ei*e 2 *) + ei*e 2 *- In 

77 



(2.7) 



and, this time, the 77-dependence of H. is explicit. This equation can be further integrated 
leading to an expression for the e-folds number N, namely 



iV-iV* 



In ~ - (1 + e u + 4* + e i*£2*) In \ + ^e u e 2 * hi 2 . (2i 



Finally, by exponentiation the above formula and by expressing the constant 0*77* in terms 
of one obtains the following equation for the scale factor itself 



a(rj) 



1 



1 + eu + e?* + ei*e 2 * - (ei* + 2e\ jf + eue 2 *) In 



1 



(2.9) 



We have reached our first goal, namely find an expression of a{rj) at second order in the 
Hubble flow parameters. 



Their derivative is indeed third order, i.e. zero at the order at which we work. 



- 6 - 



Let us now discuss how the expression of e\{ri) can be obtained. Let us notice that since 
e\ is appearing in Eq. (2.4), we need to go to third order since a term 1/ei* will remain in 
front of the final expression of Vq. This can however be obtained using the above formulas. 
Taylor expanding e\ around TV* one has 



ei = ei* + 



dei 
diV 



2diV 2 



(N - N*) 2 + 



(2.10) 



Using the fact that dei/diV = eie 2 , d 2 ei/diV 2 = e\e 2 + eie 2 e3 and the expression of the 
number of e- folds at first order [see Eq. (2.8) above], one arrives at 



ei — eu 



1 - e 2 * (1 + e u ) In (± \ + i (e 2 .* + e 2 *e 3 *) In 2 (j^j 



+ 0(e 4 ). (2.11) 



The very same method can be used to determine the second order expression of the sound 
speed. Taylor expanding c s in e-fold gives 



Zr Zr I 

c s = c s * + 



dcf_ 
diV 



(N-NJ + - 



ld 2 c 



2diV 2 



(iv - N*y + 



(2.12) 



and from the sound flow hierarchy one has dc 2 /diV = — 2c 2 5\, d 2 c 2 /diV 2 = — 2c 2 (5i(5 2 + 4c 2 <5 2 . 
Together with the expression of N — iV*, it follows that 



c 2 (??) = c s 2 



+ 2c 2 # (<5i* + (W) In ( - 4 (fc^ - 25 2 u ) In 2 ( + 0(e 3 , <5 3 ) . (2.13) 



As expected the coefficients of the logarithms are expressed in terms of the parameters 8\ 
and 5 2 . 

Finally, only the expression for 1/(77) remains to be found. By definition, one has v 2 = 
1/4 + 7] 2 U s (rj), i.e. Taylor expanding everything from the previous formulas, one gets 



V 



with 



1/(77) = v*—\ ei*e 2 * + -e 2 *e 3 * + 6uS 2 * ) In ( — ) + 0[e , 5 ) , 



_ 3 1 2 11 1 1 . 

v* = - + ei* + -£2* + 01* + ei + — ei*e 2 * + -e 2 *e3* + ei*oi* + -oi*o 2 *. 
2 2 bb 3 



(2.14) 



(2.15) 



2.3 Comoving Curvature Power Spectrum 

We have now determined explicitly the four functions appearing in the expression of the 
power spectrum V^, see Eq. (2.4). It is straightforward, although lengthy, to calculate, at 
second order, the relevant combination c 2 /(a 2 i/ei) appearing in that expression. Moreover, 
we must also find ^> . Upon using the expression of the function g(r]), one gets 



dr 



'/* 



c 2 (r)k 2 



(2.16) 
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Inserting Eqs. (2.13) and (2.14) into the previous formula and expanding everything to second 
order, the integrand in Eq. (2.16) reads 



c 2 (t)P 



V* 
T 



3 

1v* 

..2 



1 



1 . - \ 1 

£l*£2* + 2 £ 2* e 3* + "1*02* I — 



c s* 1 ( 

+ — ($1* + 61*^1*) - 1 

v* r \ 

- it {6u62 * " 26l) \ i 1 

c 4 2 1 

+ (^1* + £ i*^i*) - I 



c 2 k 2 r 2 



1 - 

1/2 



dk 2 T 2 



-1/2 



In 



cik 2 T 2 



feVln - 

-1/2 

k 2 r 2 In 2 f- 



„2 pj\ -3/2 



fe 4 r 4 ln 2 (- 

V* 



(2.17) 



Therefore, we have five different integrals to calculate in order to evaluate the term In 
the following, we write 



i=5 



(2.18) 



i=l 



and calculate each of the Ij separately. Let us also notice that the way Eq. (2.17) has been 
written is not yet fully consistent since all the terms have to be expanded to second-order. 
For instance, terms like (8i* + ei*8\*)/u* (in front of the second integral I2) should clearly 
be expanded further on in order to keep only second order expressions. For the moment, 
however, we will be keeping them this way in order to maintain clarity. Only at the end of 
the calculation these terms will be expanded. 

Let us now calculate the five integrals. Defining w = c s *kr]/iy*, which implies that 
w* = c^kij*/^* = —1, the first integral, I\, can be calculated exactly and reads 



'l-u 2 ) 1/2 +ln|d-ln 



2\l/2 



On large scales, w approaches zero and one obtains 

lim I\ = — v* (1 + In \w\ — In 2) . 

w— s-0 



(2.19) 



(2.20) 



The second integral is slightly more complicated but can also be carried out exactly. The 
result can be expressed as 



l&v* 
2 In 2 



1 . x 

£l*£2* + 2 £ 2* e 3* + 01*02h 



4 In 2 \u\ — 8 In \u\ In 



a/ 



4Li 2 



i_vT 

2 : 



u- 



(2.21) 



where Li 2 denotes the Polygamma function of order two, or dilogarithm function [66]. On 
large scales the previous expression takes the form 

_2 - 



lim I2 = 



1 



ei*e 2 * + -e 2 *e 3 * + 81*82* 4 In \w\ - 4 In 2 + — 



7r 



(2.22) 
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where we have used Li2(0) = and Li2(l/2) = 7r 2 /12 — (In 2 2)/2. We notice that 1\ and I2 
are logarithmically divergent in the limit w — > 0. We will see that this is not a problem and 
that those terms cancel out in the final expression of V^. This is expected since we know 
that the power spectrum remains constant on larges scales, and as such an exact cancellation 
of those terms constitutes a consistency check of the method. On the contrary the integrals 
^3, I4 and I5 are convergent and can be directly computed. They read 

h = v*(8u + e u 6 u )(l-hx2), h = -y {S1J2* - 25 2 u ) - 2 + 2 In 2 - In 2 2^ , 

h = j (S u + f (2 - y - 2 In 2 + 2 In 2 2) . (2.23) 

This completes our calculation of the quantity and we can now evaluate the expres- 
sion (2.4). Collecting the expressions of a, e±, c s and v established previously, one gets 
c 2 /(a 2 ve\) that has to be combined with e 2 * using the above integrals. After some lengthy 
but straightforward manipulations, one obtains 



1 + 



Hi (I8e- 3 ) 
8tt 2 Mp^uCs* 

23 4 

In2 + - In 2 2 5t + 

18 3 2 * u 



3 

2o ;2 



+ 2 In 2 ei* + 



-| + ln2 j e 2 *+ (^- ln2 



25 



7T" 

24 



^In2 + ^ln 2 2 ) 



+ 



+ 



+ 



2 ln z 2 e u 5 u + 



25 13. . 

1 In 2 

9 3 



J + jjln2-ln 2 2 

11, 1 

n 2 + - n 2 2 eL + 

18 3 2 



A3 10 



V 9 



25 



7T 



ln 2 + 2 In 2 2 ef 



1 



(2.24) 



) + 12 + i ln 2 + ln2 2 ' ei * 62 * 



1 TT 2 1 

1 h - In 2 

9 24 3 



1 



ln 2 e 2 *e 3 * 



Several remarks are in order at this stage. Firstly, in the above calculation, we have 
assumed that the initial state of the perturbations is the Bunch-Davies vacuum. This implies 
that |-Bfc| 2 = 7r/2. Notice that, in the context of k-inflation, this is a non-trivial choice 
since, as discussed in Ref. [59], the time dependence of the sound speed could be such that 
the adiabatic regime is not available anymore 4 . In this paper, we assume that this does 
not occur and that the function c s (r/) is initially smooth enough. Secondly, as announced 
above, all the time-dependent terms ln|u>| have canceled out and the expression of is 
time-independent. Thirdly, Eq. (2.24) should be compared with Eq. (51) of Ref. [59]. These 
two expressions coincide at first order, which is another indication that the above formula for 
is correct. Fourthly, in the overall amplitude, we notice the presence of the factor 18 e~ 3 . 
As explained in Refs. [48] and [59], this is typical in a approximation scheme based on the 
WKB method or its extension (such as the uniform approximation). This leads to a ~ 10% 
error in the estimation of the amplitude. In Refs. [53, 54], it was shown that, by taking into 
account higher order terms in the adiabatic expansion, this shortcomings can easily be fixed. 
In that case, one obtains a new overall coefficient which dramatically reduces the error in the 



4 Let us notice however that one can still re-define a new time variable to absorb the c s -dependence in the 
mode equation [42]. In terms of that new time variable, one could always set Bunch-Davies initial conditions 
for the scalar, but this would not be compatible with those of the tensor modes. 
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amplitude. As a consequence, we do not really need to worry about the term 18 e -3 and, for 
practical applications, one can simply renormalize it to one. 

Finally, the above expression of depends on 77* which depends on k. Our goal is now 
to make this hidden scale dependence explicit and to re-express the power spectrum at a 
unique pivot scale defined by 

k ar]o = -—. (2.25) 

This is achieved by re-writing all the quantities appearing in the power spectrum at a 
single time, r/ = rj . Technically, this means that, say, ei* should be written as ei* = 
eio — eio£2o In (??<>/??*) + C(e 3 ) and that the dependence in 770/77* should be replaced with a 
dependence in k a /k. This is performed by making use of the relation between the time 77 
and 77*: 

77* k c S o , . 

— = — v* — . (2.26) 

77o nJ C s * 

Working out the previous equation at second order, one obtains that 

. fv*\ A 3 , /, r ,2 1 2 r . r2 

In I — I = I In - + In — I II - 5 1<r - e 1<r 6u - -e 1<r e 2<r - -e 2<> e 3< , - ^io02o + o u 

2 1 2 4 1 2 4 2 1 2 

+ g£lo + g£2o + gOlo + £lo£2o — g e l« d lo + g £ 2o£3o + g<h«02o + g^U ~ Jg e 2© 

- " lt2.Su + \Sij2o In 2 I + 5uS 2 . In | In ^ + i<5 lo «5 2<> hi 2 ^. (2.27) 

This finally leads to one of the two main new results of this paper, namely the expression 
of the scalar power spectrum in k-inflation at second order in the Hubble and sound flow 
functions 



H 2 (l8e~ 3 ) [ /2 D 2 \ 

LI 1 - 2(1 + D)eu - De 2o + (2 + D)S U + [g+ D + ^) S l 



%TT 2 Mleuc, 



+ + 2D + ^ - ^ ^^o + (~ -3D- 2D 2 ) e lo 5 u + (y + 2D + 2D 2 , , lo 

+ (I " D ~ D2 ) ^ + ("T " ^ + ^ + s) — + (I + T) * 

+ ( 24 " 18 ~ ~2~ ) €2<>e3 ° + [~ 2€u ~ £2 ° + 6u + (1 + D)<y2 ° 

+ (2 + £)<5 lo <5 2o " (3 + 4D)eu5u + 2(1 + 2D)ej - (1 + 2 J D)e 2 o<5i - (1 - 2D)e u e 2<> 
k 

ka 



k ( 11 11 

+ De 2<> - £>e2o£3o] hi — + 2e 2 u + eue 2<> + -e| - -e 2 o£3o + r#i, + -^lofo 



-2e 1 o ( 5 1 o-e 2 o5io) ln 2 ^|, (2.28) 

where we have introduced the quantity D defined by D = 1/3 — In 3. One easily checks that, 
at first order, this expression exactly coincides with Eq. (53) of Ref. [59]. More details in the 
comparison of the above formula with the existing literature can be found in Sec. 3. 
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2.4 Tensor Power Spectrum 

In this section, we repeat the previous analysis but for tensor perturbations. Since the 
method is the same and, fortunately, the calculations are easier, the details will be skipped. 
The main difference between gravity waves and density perturbations is that their effective 
potential is not the same, see Eqs. (1.4) and (1.6). This implies that the function 1/(77) f° r 
tensors is different from the one of the scalars. One gets for the tensor 



v (V) = 1 + 3ei* + 4ef, + 4ei*e 2 * - 3ei*e 2 * In 1^-1 + 0{e 6 ) . (2.29) 

As a consequence, the functions g(rj), f(rj), and hence ^f, are also different. Using the uniform 
approximation to evaluate /i^ and inserting the corresponding formula into the expression of 
Vh given by Eq. (1.9), one obtains 



_ 2(l8e- 3 )i^ 2 

h 7T 2 M 2 



-^ + 21n2lei*+ I ^— - y + | In 2 - In 2 2 ) 



,13 10, , o 

(___l n2 + 21n 2 2 ) 6 2 , 



(2.30) 



This equation is for the tensors what Eq. (2.24) is for the scalars. As explained before, one 
has still to make explicit the scale dependence hidden in 77*. In the case of tensors, the pivot 
point is usually defined by k+rj+ = — 1 since gravity waves propagate at the speed of light. 
This leads to the following expression for the power spectrum 

2(l8e- 3 )# 2 f , , /17 2 \ 2 / 19 7T 2 

V h = \ 2M l - {l - 2(1 + D)^ + (j + 2D + 2L 2 ) 4 + (- T + ^ 

-2D- D 2 ^j ei*e 2 * + [-2e lir + 2(1 + 2£>)ef* - 2(1 + D)e u e 2 *] In ^ 

+ (2eL- ei , e2 ,)ln 2 ^|. (2.31) 

Let us notice that, in order to obtain this relationship, we have used the initial conditions for 
gravity waves \Bk\ = 1/M 2 j. Otherwise, one notices the presence of the WKB factor 18 e 3 
and one can check that, at first order, it coincides with the known expression for the tensor 
power spectrum. The above formula, being expressed at the time 77*, is convenient for SFMK 
models only, but not for k-inflation. Indeed, all parameters here are functions evaluated at 
the time 77* which is different that the one at which the scalar power spectrum is calculated, 
namely t? . It has become a common mistake to try fitting data with both Eq. (2.28) and 
Eq.(2.31) while implicitly assuming that all Hubble and sound flow "parameters" are the 
same. As we have explicitly shown before, they do differ and such a fit would absolutely 
make no sense. 

However, within slow-roll, one can re-express the tensor power spectrum at the same 
pivot point as for the scalar power spectrum. As before, each quantity in the tensor power 
spectrum should be re-expressed at the scalar pivot point, as for instance ei* = ei — 
ei e 2o lnc s<> + 0(e 3 ,<5 3 ). The quantity c so appears because it is present in the ratio of the 
tensor to scalar pivot points. It follows that the final expression for the tensor power spec- 
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trum for k-inflation is 



2 (I8e~ 3 ) Hi . 
Vh = ^ { 1 - 2(1 + £> - lnc so )e lo + 



7T 2 M 2 



17 



+ 2D + 2D 2 + 2 In 2 c s 



2{l + 2D)lnc s 



4 + 



19 7T 



+ — - 2D - D 1 + 2(1 + D) lnc s<> - In 2 c s 



+ 



9 12 

-2e u + {2 + 4D- 41nc so )eL + (-2 - 2D + 2 In c so )e u e 2c 



+ (2e? - e lo e 2o ) In 2 



hi' 



eio£2o 



(2.32) 



where now "diamonded" terms are evaluated at the scalar pivot point. This new formula 
is the second main result of the present paper. It extends to second order the results of 
Ref. [59]. 



3 Discussion and Conclusions 

The power spectra of Eqs. (2.28) and (2.32) represent the main result of this article. There 
are the first calculation, at second order in the Hubble and sound flow functions, of the scalar 
and tensor power spectra in k-inflation within the uniform approximation. In this section, we 
discuss our results and check their consistency. In particular, in some limits, our calculation 
should reproduce known results already derived in the literature. As we show below, this is 
indeed the case. 

We have seen before that the power spectrum is obtained as an expansion around the 
pivot scale and that the most general expression of can be written as 



(3.1) 



n>0 



where is the overall amplitude and the coefficients a n are functions of the horizon flow 
parameters. The expression of a n always starts at order n, i.e. ao starts with one, a\ 
starts with a term of order 0(e, 5), a 2 with a term of order C(e 2 , <5 2 , e<5) and so on. As 
already mentioned before, k-inflationary power spectra were determined at first order in 
Ref. [59]. This means that the expression found in that paper included only the first two 
terms, proportional to ao and a%. There is however a trick derived in Ref. [67] which allows 
us to determine some higher order terms. Indeed, the power spectrum should not depend 
on the choice of the pivot scale, which is arbitrary. As a consequence, one can establish the 
following recursion relation 



dlnP c 
din k a 



a n + 



da r 



din k a 



(3.2) 



Given that oo was given at first order in Ref. [59], it was then possible to calculate a\ up to 
second order and a 2 to third order [see Eqs. (64) and (65) in that reference]. Therefore, one 
can compare those formulas to the expression obtained in this article. One finds that they 
are the same, indicating the consistency of our results. 

Another way to verify the validity of our expressions is to take the limit c s = 1 and to 
compare the resulting formulas to the results already obtained in the literature for SFMK 
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models. As mentioned in the introduction, second order results were first obtained using 
the Green function method in Ref. [47]. The corresponding expression for the scalar power 
spectrum reads 

(n2 „ 77T 2 \ fC 2 7T 2 \ 2 / C 2 7T 2 \ 

+ { c - c + l2- 7 )^ + {-2 + Y- 1 ) e2+ {-T + u)^ 

+ [-2 £l - e 2 + 2(2C + l)e\ + (2C - l)e 1 e 2 + Cei - Ce 2 e 3 ] In A. 
+ ( 2e\ + e 2 e 2 + \e\ - \e 2 e 3 \ In 2 A I, (3.3) 



2 z 2 "/ fc, 



where the constant C is defined by C = 7 + In 2 — 2 ~ —0.7296, 7 being the Euler constant, 
while the expression of the gravity wave power spectrum can be written as 

n = - 2«7 + Der + (2C 2 + 2C + £ - 5 ) + (-C 2 - 2C + g - 2) £l£2 

+ [-2 £l + 2(2C + l)e 2 - 2(C + l)eie 2 ] In + ^2ej - e x e 2 ^j In 2 |- 1 . (3.4) 

In the two previous formulas (3.3) and (3.4), the Hubble flow functions are evaluated at 
time 77® such that a(r] m )H(r] m ) = k m which slightly differs from the time k r] = — 1 (for 
c so = 1) used in the present paper. Therefore, if we want to compare Eqs. (2.28) and (2.32) 
with Cs = 1 to Eqs. (3.3) and (3.4), one should first re-express the latter in terms of the 
Hubble flow parameters evaluated at time k^^ = — 1. In the following, in order to simplify 
the discussion, we focus only on the scalar case but the tensor case could be treated in the 
same manner. From the definition of rj m one has r] s> /r] = 1 + ei + ef + ei e 2<> + 0(e 3 ). As 
consequence, in Eqs. (3.3) and (3.4), one should just replace e-y, e 2 with ei , e 2o and H 2 je\ 
with H 2 /eu(l + 2e 2 + eue 2o ). This yields the following expression 

m f .. ( -, 2 ... - 2 



= ^ 2 m ~ 1 1 - 2(C + l)6i. - Ce 2o + \2C 2 + 2C + - - 3 j ef 
, 7vr 2 e \ , , /C 2 , vr 2 \ , / C 2 vr 



+ (C* - C + — - 6) e u su + ^— + — -lje z 2<> + \- — + —j e 2 ^ 
+ [-2e u ~ £2. + 2(2C + l)e 2 + (2C - l)ei e 2o + Ce 2 ^ - Ce 2a e 3o ] In A 
+ ^2e 2 + ei e 2o + ^e 2 , - In 2 A | , (3.5) 

that can be now compared to Eq. (2.28). As already discussed, the overall amplitude differs by 
the WKB factor 18 e~ 3 . We also notice that the terms in D in Eq. (2.28) exactly corresponds 
to the term in C in Eq. (3.5). For instance, the coefficient of ef in Eq. (2.28) contains a term 
2D 2 + 2D while the coefficient of e 2 in Eq. (3.5) contains a 2C 2 + 2C. One easily checks that 
this is the rule for all first and second order terms. Provided one substitutes D with C, the 
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first order term in the amplitude, the coefficient of lnfc/fc and the coefficient of ln 2 (fc//c ) 
are identical. The only difference appears in the second order terms in the amplitude. For 
instance, the coefficients of ef in Eq. (2.28) is 2D 2 + 2D + 17/9 while it is 2C 2 + 2C + vr 2 /2-3 
in Eq. (3.5). But 17/9 ~ 1.88 and 7r 2 /2 — 3 ~ 1.93 and, therefore, the two terms are in fact 
numerically very close. The same is true for all the other terms in the amplitude. Therefore, 
we conclude that our result (2.28), specialized to SFMK models, is fully consistent with 
Eq. (3.5) that comes from another approximation scheme. This confirms its validity. 

Let us now compare our result to the one of Refs. [53, 54] calculated with the help of 
the WKB approximation. The scalar power spectrum obtained in those articles reads 

H 2 , , „ f o _ l\o 



n = 87r 2 M 2 / wkb 1 - 2(Avkb + " D WKB e 2 + I 2D^ + 2D 



- ei 



7T 2 20 \ (D* 2\ o f TT 2 



+ { D ^- D ^ + T2-y)^ + {^ + 9)^ + [-^ KB 
+ [-2 ei - e 2 + 2(2D WKB + l)e 2 + (2Z? WKB - l) Cl e 2 + D WKB e 2 - D WKB e 2 e 3 ] In 

+ (2e 2 + e l£2 + \e\ - ±e 2 e 3 ) In 2 ^ | , (3.6) 

while the tensor power spectrum is given by the following formula 

Vk = 1 1 " 2(j ° WKB + 1)61 + (^WKB + 2 ^WKB " l) 4 + ("^WKB " 2 ^WKB 




7T 2 19\ fc 

12 " Y J 6162 + [ _2ei + 2 ( 2D wkb + " 2 Pwkb + l)eie 2 ] In ^ 



(2 £ 2 - eie2 )ln 2 ±i. (3.7) 



h 

n>® J 

In these equations, A WKB = 18 e -3 and D WKB = 1/3 — In 3, that is to say exactly what 
was found by means of the uniform approximation as D WKB = D. As already mentioned, 
Refs. [53, 54] have shown that, by taking the next order in the adiabatic approximation 
into account, one obtains a new value for these two constants (in some sense, they are 
renormalized) , namely A WKB becomes 361/(18e 3 ) ~ 0.99 and D WKB = 7/19 — In 3 ~ —0.7302. 
In particular, the new value of D WKB is closer to the constant C than the non-renormalized 
one. Both Eqs. (3.6) and (3.7) are evaluated at the pivot time ry® and have to be time-shifted 
to r] to be compared with our results. Proceeding as previously, it is easy to show that this 
modifies the coefficients of e 2 which now becomes 21?^^ + 2D WKB + 17/9, and the coefficient 
of eie 2 which becomes 2Z) 2 VKB — D WKB -j-7r 2 /12 — 11/9. In other words Eqs. (2.28) and Eq. (3.6) 
expressed at rj are exactly the same for c so = 1. This is maybe not so surprising considering 
the fact that the WKB and uniform approximations are closely related methods. 

A few words are in order about Ref. [55]. Historically, this is probably the first paper that 
attempted to evaluate the k-inflationary power spectrum at second order in some equivalent 
of the Hubble and sound flow functions used here. The method chosen is the Green function 
expansion discussed before. However, a specific form for the sound speed, which in the 
language of the present paper would be a first order approximation of c 2 , was also postulated. 
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Together with a k— dependence kept implicit, this makes the comparison with the present 
work difficult. For this reason, we do not investigate further this issue. 

To conclude, let us briefly recap our main result and discuss directions for future works. 
In this paper, using the uniform approximation, we have calculated the scalar and tensor 
power spectra in k-inflation, at second order in the Hubble and sound flow parameters, see 
Eqs. (2.28) and (2.32). We have carefully checked that, in the various limits where our 
calculation reduces to known cases, consistent results are obtained. The next step is clearly 
to use these power spectra in order to constrain the values of the Hubble and sound flow 
parameters using CMB observations. This was done in Ref. [62] but only for the first order 
power spectra (since only this result was available at that time). Given the on-going flux 
of high precision data, such as those from the Planck satellite, the results obtained in this 
article should be important to keep theoretical uncertainties at a minimal level. In this way, 
as discussed in the introduction, one may hope to obtain unprecedented information on the 
inflationary scenario. 
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